Entanglement and its evolution following a quench in the presence of an energy current 
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We study the Ising spin chain with a Dzyaloshinskii-Moriya interaction focusing on the static and dynamic 
properties of the entanglement entropy, following a quantum quench. We show that the effects of the additional 
anisotropic interaction on the phase diagram and on the dynamics of the system are captured by the properties 
of the entanglement entropy. In particular, the model provides a way to study the quench dynamics in a system 
with an energy current. We consider quenches starting from an initial excited state of the Ising spin chain, and 
we analyze the effects of different initial conditions. 



I. INTRODUCTION 

Recent experiments have shown that it is possible to study 
unitary nonequilibrium dynamics in quantum systems on long 
time scales [1]. One important result of these investigations 
has been the observed lack of thermalization. The reason for 
this non-thermal behaviour is attributed to the near integrabil- 
ity of the system. This observation has motivated recent stud- 
ies of the non-equilibrium properties of integrable quantum 
model Hamiltonians [2]. The interest in this topic is also mo- 
tivated by its relevance to a variety of experimental situations, 
including cold atoms [3], Penning traps [4] and Josephson- 
junction arrays JD]. 

There are many ways in which it is possible to drive a sys- 
tem out of equilibrium. Quantum quenches, and coupling to 
baths with different temperatures (or potentials) are the most 
common protocols. While in the quench scenario the focus 
is on the unitary dynamics of the full system, in the second 
case one usually deals with an effective description of the dy- 
namics of the subsystem only. Different out-of-equilibrium 
dynamics have nevertheless similar characteristics, for exam- 
ple the presence of currents (of particles, energy, or heat) . 

In this work, we study an Ising-like model system driven 
out of equilibrium with a quantum quench in the presence of 
an energy current. In the usual quench protocol the system is 
prepared in the ground-state in the absence of any current, and 
subsequently the dynamics drives it to some excited state. The 
model Hamiltonian we consider here allows for the study of 
two different new situations: a quench from the ground-state 
in the presence of an energy current, and a quench from an ini- 
tial excited state in the absence of an energy current. The two 
scenarios are associated with different Hamiltonians which, in 
a sense, are dual to each other. To characterize the system's 
behaviour we focus our attention on the Entanglement En- 
tropy (EE), as measured by the Von Neumann entropy, which 
is a central quantity in the characterization of nonequilibrium 
quantum dynamics Q6[|. 

The structure of the paper is the following: in Sec. II we 
introduce and describe the properties of the model Hamilto- 
nian, in Sec.III we consider the static properties of EE for this 
model, and in Sec. IV we study the dynamics of entanglement 



following a quench. Finally, we states our conclusions. In the 
Appendix the reader can find more details of the calculations. 



II. THE MODEL 

We consider an Ising spin chain in transverse field Hi, with 
an additional Dzyaloshinskii-Moriya (DM) interaction Hdm- 
The total Hamiltonian Hj + Hd m is defined as follows 
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where h is the external magnetic field, and ( is the cou- 
pling parameter determining the strength of the DM interac- 
tion. Such an anisotropic interaction is present in many low- 
dimensional materials with the necessary crystal symmetry, 
and it originates from spin-orbit coupling |]7l|8|]. Furthermore, 
the DM interaction is of relevance in quantum information 
theory, since it plays an important role in the physics of quan- 
tum dots H, and in fault-tolerant quantum computation II 1 CHI - 
Adding the DM term to Hi does not affect the solvability 
of the model Hill , and interestingly enough it provides the 
system with a richer phase diagram. These features have been 
used in [ 12] to study the effective out-of-equilibrium quantum 
dynamics of the model. Hdm can be viewed as a current 
term. The reason for this is the following. The equation of 
motion for the local energy density of Hi, defined by ej — 
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One can write the time derivative of the energy current as the 
divergence of the energy current 
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It thus follows that Hdm is precisely the sum over all sites of 



the local currents ^ . Cj . 
tation value of Hdm 



Therefore, the ground-state expec- 
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becomes an order parameter indicating the presence of an en- 
ergy current. Once the total Hamiltonian has been diagonal- 
ized, which can be done with the usual Jordan- Wigner and 
Bogoliubov transformations 

H = Y / ^ q b%, (6) 

the effect of the DM interaction is clearly observed at the 
single-particle level. The single-particle spectrum is given by 

A q = \f\ + h 2 + 2h cos q + ( sin q, (7) 

with q G [— 7T, 7r) the momentum of the quasi-particle. As 
can be seen in the above expression and in Fig[T] the DM in- 
teraction makes the spectrum non-symmetric with respect to 
q = 0. Note that the ground state of the Hamiltonian includ- 
ing the DM interaction is the same for all values of £ in the 
interval [0, 1]. In particular this means that, within this range 
of values, the Ising model in a transverse field Hj and the sys- 
tem described by EqQ]have the same ground state. Beyond 
C = 1 (with h < the Fermi sea starts to be populated by 
the modes in between the zeros of the single particle spec- 
trum (i.e. between q + and q_ in FigJTJ. This implies that the 
ground state is not anymore that of Hj. Furthermore, since 
the DM term commutes with the rest of the Hamiltonian, at 
the many-body level when £ = 1 we must have a level cross- 
ing between the ground state of Hi and some previously ex- 
cited Hamiltonian eigenstates. Fig|2]shows the phase diagram 
of the model 11211 . There are three regions: ferromagnetic, 
polarized paramagnetic, and the so-called current phase, char- 
acterized by J 7^ 0. The current phase is gapless, and the 
two-point correlation functions show a power-law behaviour 
with an oscillatory amplitude: {&f&f+ n )gs ~ cos(fcn), 
where Q is a non-universal function and k = arccos | lfl2tl . 

In the following, we study the entanglement properties of 
this model, and subsequently we analyze new quench proto- 
cols for this spin system. 

HI. STATIC ENTANGLEMENT ENTROPY AND PHASE 
DIAGRAM 

In this section we show how entanglement can be used to 
characterize the different phases of the model. To measure 
the EE, we consider a bipartition of the spin chain into two 
subsystems A and B. For this setup a good measure of EE be- 
tween the two partitions is given by the von Neumann entropy 
Sa = — TrpLln pl, where p^ is the ground-state reduced 
density matrix of the subsystem A with L spins. 

It is known that for critical one-dimensional systems the EE 
scales logarithmically in the subsystem size, with a prefactor 
given by the central charge of the associated Conformal Field 
Theory (CFT), 

S L = ~\nL + S , (8) 

where c is the central charge and So is a non-universal con- 
stant lfl3l Q • On the other hand, in the non-critical region 
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FIG. 1. (Color online) Spectrum of the Hamiltonian in Eq[T] We 
show the spectrum for 4 different values of f, while keeping h = 0.5 
fixed. See also figures in 11211 . 
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FIG. 2. (Color online) Phase diagram in the h — £ plane of the model 
in Eq[T] The dotted line is a critical line where the model shows the 
same universal properties of the quantum Ising model in transverse 
field. See also figures in lfl2ll . 



of the phase diagram, the entanglement entropy saturates to a 
value which depends on the correlation length £, 

5t a | Inf. (9) 

Both Eqj8]and Eq|9] characterize the ground-state properties 
of EE for one-dimensional systems. 

Apart from the ground state it is also of interest to investi- 
gate entanglement properties of excited states. Recently, two 
works have appeared on this topic. In lfl5ll it has been shown 
that there are excited states for which the logarithmic scaling 
of EE can have prefactors different from the ground state, and 
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FIG. 3. (Color online) Sl vs L at h — 1.0 for different £. The scaling 
behavior changes from Sl ~ | In L on the critical line separating 
the ordered and disordered phase to Sl ~ 3 In £ inside the current- 
carrying phase. 



that for some excited states the scaling can be extensive in 
the subsystem size, instead of logarithmic. In [16] the authors 
have studied the connection between EE for excited states and 
properties of the associated CFT not contained in the central 
charge. The EE of excited states is of interest also in the 
context of quantum quenches since in this setting the system 
is unitarily driven from the initial ground state to an excited 
state. 

The Hamiltonian in Eqfl] naturally fits in this set of prob- 
lems. Following the discussion from the previous section the 
model we are considering allows us to study the EE of some 
excited states of Hj, simply by tuning the coupling constant 
associated with the DM term. The excitations we can consider 
in this way are characterized by the modes in between the ze- 
ros of the single particle spectrum that get populated when 
C > 1 and h < ( ©. For these states the EE can be evalu- 
ated analogously to the ground state of Hj (see also B17I1 for 
related analytical study). 

Let us consider in detail the entanglement properties of the 
different phases shown in Figj2] First, we compare the scal- 
ing of the entanglement in the non-current-carrying critical 
regions and in the region where an energy current is present. 
Fig|3] shows the result of the simulations for the scaling be- 
haviour of the ground-state EE with £ < 1 and h = 1. For all 
values of £ 6 [0, 1] one observes the same scaling result. In 
the non-current-carrying region critical states are present only 
on the h = 1 line of the phase diagram. On this line, sepa- 
rating the ferromagnetic and polarized paramagnetic phases, 
the ground state of the system is the same as in Hi. This im- 
plies that the EE scaling is logarithmic with a prefactor of c/3, 
and c = 1/2. Note that also the entire current-carrying phase 
(C > 1 and h < Q is gapless, and in this sense critical. At any 
point in this phase we observe logarithmic scaling of the EE 
in the subsystem size. This is consistent with the discussion 
in the previous section on the algebraic decay of the two point 
correlation function. Interestingly, the prefactor of the loga- 
rithmic scaling of EE in the current phase is twice as large as 
the prefactor in the non-current phase. This doubling reflects 



the increased number of zeros in the single-particle spectrum, 
consistent with the results of [17[]. In fact, when £ > 1 the 
ground state of EqQ] is the filled Fermi sea of modes in be- 
tween <7_ and q + (see FigQ}. Since the ground-state in the 
current phase is effectively an excited eigenstate of Hj, we 
could expect a scaling of EE that is extensive in the system 
size, as shown in Il5ll for excited states. The reason why this 
is not the case is due to the nature of the single-particle spec- 
trum, which at most can have two zeros (see FigJTJ, and thus 
does not satisfy the requirements found in IU5I1 for an exten- 
sive scaling of EE. 

The DM interaction in the Hamiltonian affects also the sub- 
leading term in the scaling of EE Sl = 5 InL + So(h, £). 
Deriving the analytical form of the sub-leading order term 
So(h, C) is complicated. Nonetheless, one can investigate this 
term numerically. In Figj4]we see that Sq is constant on the 
critical line h = 1 with ( < 1. As soon as ( > 1 and h < (, 
So increases, but becomes almost constant for large £. Also 
So is maximum at h = 0, and So is minimum at the critical 
point h = £. From the behaviour of So we can conclude that 
a given block has the highest entanglement when all the neg- 
ative modes {q 6 [— 71", 0)) in the Fermi sea are filled. Conse- 
quently EE increases with higher values of the energy current. 

We now consider the differences between the critical lines 
shown in the phase diagram (Figj2j, separating different 
phases. The only second-order quantum phase transition is 
found along at the h = 1 line (with ( < 1), which corresponds 
to the Ising quantum phase transition (see Fig. |5}. On the 
hand the boundaries of the current-carrying phase with both 
the paramagnetic and the ferromagnetic phases are character- 
ized by a level crossing. This translates into a sudden jump in 
EE (see Fig|5]and Fig|6]l. The value of EE is always higher 
in the current carrying phase because of the presence of long- 
range correlations that decay algebraically. The plots in FigO 
and Fig|6]show that controlling the DM term can be used as 
an entanglement switch. The amount of entanglement can be 
driven by the DM coupling term or the magnetic field, which 
are controllable parameters in optical lattices 111 811 . 



IV. ENTANGLEMENT DYNAMICS FOLLOWING A 
QUENCH 

In this section, we focus on the quench dynamics of the 
EE. Quenching provides a way to excite a system, initially 
prepared in the ground state, and to subsequently study the 
non-equilibrium dynamics of the model (in the following, we 
denote with a subscript the value of the parameters describ- 
ing the initial Hamiltonian). As stated previously, the model 
in EqQ]is of interest because it combines two different mech- 
anisms typically used to drive a system out of equilibrium: 
quantum quenching, and the coupling to a field originating 
a current in the system. Furthermore the inclusion of the DM 
term allows us to study a model Hamiltonian where the energy 
current can be controlled and used in the quench protocol. 

In our setup, the quench can either involve the magnetic 
field h, the DM coupling ( or a combination of the two. Since 
the DM term commutes with the Hamiltonian, a quench in ( 
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FIG. 4. (Color online) Non-universal nature of So. (upper panel) So 
vs. £ at different h; (lower panel) So vs. h at different £ . 




FIG. 5. (Color online) Sl vs. h at L = 60 for different £. (up- 
per panel) Static entanglement along the transition between the or- 
dered ferromagnetic and the disordered paramagnetic phase. The 
peak signals the presence of long-range correlations at the critical 
point, which is a signature of a second-order quantum phase transi- 
tion, (lower panel) Static entanglement along the transition between 
the disordered paramagnetic and the current-carrying phase. The 
sudden change in entanglement followed by the absence of a peak 
at the critical point is a signature of the first-order quantum phase 
transition. 



leaves the system in one of its eigenstates, providing a trivial 
evolution of the EE. On the other hand, quenches in the mag- 
netic field give more interesting behaviours. If the quench is 
done with the initial state prepared in a region with no current 
the results are similar to those found in [6], where quenches 
for the Hj were considered. This is due to the fact that, in 
the absence of a current, the ground-state wave function ini- 
tially is identical to that of Hj, and the time evolution is not 



FIG. 6. (Color online) Sl vs. ( at L = 60 for different h. The plots 
for h — 0.5 corresponds to the static entanglement along the tran- 
sition between the ordered ferromagnetic and the current-carrying 
phase. This is a first order quantum phase transition which occurs at 
£ = 1 (analogously for h = 1). The plots for h = 2 and h = 3 cor- 
respond to the static entanglement along the transition between the 
disordered and the current-carrying phase. This is a first order QPT 
which occurs at f = h. 



affected by the presence of the DM term (see the derivation of 
Eqf25]in the Appendix for a proof of this). More interestingly, 
if the quench involves an initial state inside the current phase, 
new non-trivial behaviours can be expected, since the ground 
state now is radically different. It is important to notice that 
the DM coupling enters only in the specification of the ini- 
tial state, whereas the evolution can be effectively described 
by the Hamiltonian without the DM term. The calculations 
showing that this is in fact the case can be found in the Ap- 
pendix. 

We first compare the evolution of the EE for different 
quenches inside the current-carrying phase. Fixing the cou- 
pling constant of the DM term, and quenching only the exter- 
nal magnetic field we obtain the results shown in FigJT] One 
always has an initial ballistic evolution of the EE, which grows 
linearly in time (measured in units where the speed of the ele- 
mentary excitation is unity) and saturates at some point. Quite 
interestingly, the saturation time (hence also the rate at which 
entanglement is initially building up) depends on the particu- 
lar evolving Hamiltonian. This way we can control the time 
needed to generate the maximal asymptotic amount of entan- 
glement. This property is relevant also from a computational 
point of view. In fact, DMRG-like schemes, used for the sim- 
ulation of the time evolution of quantum systems, can take 
advantage of the lower rate at which entanglement is gener- 
ated. Knowing the regions in the phase diagram where such 
rates are lower can provide more efficient time simulations. 
As far as we know this is a new feature that is not present in 
other quench protocols considered so far in the literature. The 
other aspect that is important to notice in FigJ7]is the special 
role played by the line h — 1 in the phase diagram, which 
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FIG. 7. (Color online) Quenches from the current carrying phase. 
S L (t) vs. the time steps, with L = 60, h = 4.0, Co = ( = 5.0 for 
different h. Note that the extent of the initial linear regime depends 
on the particular evolving Hamiltonian. 



turns out to provide the maximum asymptotic EE for different 
quench parameters. This can be understood by mapping the 
quench for Hj + Hdm to a quench protocol for Hi only. As 
stated in the previous section, the entanglement evolution with 
respect to H(h, () is identical to the evolution with respect 
to H(h,0). Furthermore, the ground state of H(Hq, Co), the 
initial Hamiltonian in the quench protocol is also an excited 
eigenstate of H(ho,0), because of the commutativity of the 
DM term with the total Hamiltonian. From this dual perspec- 
tive the effect of the current is that of effectively quenching 
an excited eigenstate without the current term. For the Ising 
model, a quench from ho ^ 1 yields the maximum value of 
Sl(oo) when quenched to h = 1, because the energy gap 
closes at h = 1, and hence a large number of zero energy ex- 
citations can be produced. While the asymptotic value of EE 
depends on the particular excited state at the beginning of the 
quench. 

Fig|8] shows results of simulations for quenches with in- 
creasing values of the DM field in the current-carrying phase. 
The asymptotic value of the EE decreases with increasing (. 
This is consistent with the phenomenological picture provided 
in |@], and with the fact that if the system starts in an excited 
state, the available number of unoccupied modes that can be 
occupied after the quench is smaller than in the case of having 
the ground state as an initial state. Furthermore, Figj8]shows 
that the time at which the EE saturates does not depend on 
C, and consequently does not depend on the particular initial 
Hamiltonian eigenstates (as long as it is not an eigenstate of 
the evolving Hamiltonian). The line with £ = 100 in FigJH] 
shows that very deep into the current phase quenching does 
not create entanglement. In fact, when £ ^ ho quenching the 
magnetic field is just a small perturbation to the Hamiltonian, 
which then approximately stays in the ground state. 

Finally, we verify that the presence of an energy current 
does not affect the extensive nature of the asymptotic value 
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FIG. 8. (Color online) Quenches from the current carrying-phase 
with different values of the current driving field C- vs. the 

time steps with L — 60, ho = 3.0 to h = 1 for different Co = C- 
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FIG. 9. (Color online) Sx(f) vs. the time steps inside the current- 
carrying phase for different block sizes L. Quenching is done from 
h = 2.0 to h = 1.0 with Co = C = 3.0. 



of EE (Fig|9]l, and its proportionality with the quench size 
(FigO. 



V. CONCLUSIONS 

We have studied the static and dynamic properties of the 
entanglement entropy in the Ising spin chain with a transverse 
field and a Dzyaloshinskii-Moriya interaction. The model is 
characterized by the presence of an energy current for certain 
regions of the phase diagram. 

Concerning the static properties we have analyzed the tran- 
sitions between phases with no energy current and the phase 
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FIG. 10. (Color online) Seo(t) vs. the time steps for quenches to h = 
1 from various ho inside the current-carrying phase at £ = £o = 4. 

where an energy current is present. The transition is captured 
by a discontinuity of EE as a function of the parameters, and 
by a distinguishable scaling behaviour in the current-carrying 
and non-current-carrying regions. In particular, the leading 
logarithmic term of the EE scaling with respect to the system 
size has a prefactor in the current-carrying region which is 
twice as large compared to the second order Ising critical line. 

Concerning the behaviour of the entanglement evolution 
following a quench, the model in Eq [Hallows us to study new 



quench protocols. The usual schemes consider quenches from 
an initial ground state. This scenario, for the model in Eq. fl] 
effectively corresponds to a quench from an initially excited 
state of the Ising spin chain in transverse field (without DM 
interaction). The main result of this analysis shows that the 
ballistic picture presented in is still valid, although with a 
significantly different aspect. In particular the entanglement 
saturation time in the current-carrying phase depends on the 
details of the evolving Hamiltonian. This is an indication of 
the role played by the evolving Hamiltonian on the propaga- 
tion of excitations. This result is of relevance in tuning the 
dynamics of the system in regions with a different rate for the 
propagation of entanglement. Furthermore it also provides a 
characterization of the regions in the phase diagram that can 
be simulated more efficiently with DMRG-like techniques. 

From a general point of view, the model in Eq. [T]also sug- 
gests a simple way to study the quench dynamics of initial 
excited states in integrable systems. The addition of a com- 
muting term in the Hamiltonian causes a reshuffling of the 
spectrum that, without changing the integrability of the model, 
allows us to obtain non-trivial results about the excitations in 
the original model. The same trick can be applied to other 
systems of interest. 

We thank Letian Ding, Zoltan Racz, and Paolo Zanardi 
for useful comments. This work has been supported by NSF 
grants PHY-803304, and DMR-08049 14. 
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Appendix. In this Appendix we give a detailed descrip- 
tion of the steps involved in first evaluating the entangle- 
ment entropy of the Hamiltonian in EqQ] and then calcu- 
lating its time evolution. After the standard sequence of 
Jordan- Wigner and Bogoliubov transformations the Hamil- 
tonian is in the diagonal form H = YTk=--K ^ko\,bk, with 
A k = + h 2 + 2/icosfc + Csinfc). 

The density matrix of the subsystem od size L, embedded in 
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a system of size N, can be obtained tracing out the rest of the 
system 



PL = Tr N - L (p) = A e u , 



(10) 



where Aq is a normalization constant and H is a quadratic 
hermitian operator 

L 1 

"H can be diagonalized via a generalized Bogoliubov transfor- 
mation. The reduced density matrix has the form 



p L = A exp[- ^2 £ q d l d q\ ( 12 ) 

9=1 

UsingTr( j o L ) = 1, weget A = U^ =1 1+c ^ ( _ £q) . This gives 
the final form of the density matrix as 



pl = nfu 



L exp(-£gdld q ) 
q=1 1 + exp(-e,) ' 



(13) 



l-cxp(-e g ) 
l+cxp( — e q ) ' 



Defining i/ 9 = 1 — w e can write 



_ / -±-2- 



(14) 



and also 



PL = n^ =1 i±-l exp[- m (i±^)4d g ] 
Z 1 v n 



1 + ZA, 



L 
9=1 



(15) 



Using the fact that cxp(d} q d q In A) = 1 + (1 — X)d q d q , one has 
for the entanglement entropy 

Sl = -Tr(p L \n(p L )) 



^[ln(l+exp(-e 9 )) + T 



9=1 



+ exp(e g ) 



^(i-^lni±^ + i-^ln---i). (16) 



9=1 



We have to calculate v q , from which we can obtain the block 
entropy. v q is given by the expectation value of d q d q and d q d q 

exp(-e„) 1 - v q 



(dgdq) = Tr(p L d\dq) = - 
{d q d\) = Tr{ PL d q d\) = — — L 



exp(-e g ) 2 

l + Va 



q ql KI L q q > l + exp(-e g ) 2 



(17) 



We define four 2L x 1 column vector: 7J 



, D 



- / V ,/ ) a,Kl I c 



it 



where 



i = (ii, . . . , d£) , and similarly for c. The previous Bogoli- 
ubov transformations can be expressed in a compact matrix 
notation as 



D=\{ h g \C. 



and 



D = C* 



9' 
h' 



9 



(18) 



(19) 



where g and h are LxL matrices. In terms of expectation 
values we have 



(DD 



g h 
h g 



9" h 



(20) 



Let us now consider a quantum quench protocol. Initially 
the system is prepared in the ground state of an Hamiltonian 
W , and suddenly one of the parameters is changed, and the 
new Hamiltonian is denoted by H . The quasi-particle operator 
vector B' k = (b' k , b 1 _ k ) 1 i s associated with H' , and the vector 

Bk = (bk, b_ k Y is associated with H. Similarly for Ck, 
which is associated with bare vacuum fermions. Define also 
the matrix i? M (a) = cos(§)I + ia^ sin(^), where <r M are the 
Pauli matrices, and p = x,y, z. It can be easily seen that 
Cfc = R x (dk)Bk, and Ck = R x (8 'u}B k , with Bk a parameter 
of the Bogoliubov transformation lfl9ll . From which we can 
write B k = R x (9' k - e k )B' k . 

When a quench takes place, the time evolution of the 
fermion operators is given by Bk(t) = e~ %Ht Bke . We 
can write Bk(t) = S z (—2Akt)Bk, where 



S z (-2A k t) 



zAibt 









-iA- h t 



Notice that the energy spectrum of the Hamiltonian in Eq{T]is 
not symmetric, which means that in general A_fc ^ In 
order to evaluate the two-point correlation functions we con- 
sider different cases. When the initial state of the system is in 
the non-current-carrying region we have 



(B k B' k ) 



1 





(21) 



If the initial state of the system is in the current-carrying phase 
and k g (kx, kz), where fci and &2 are the zeros of the spec- 
trum then 








For k lying between —k\ and — k-2 we have 



(B'kB'k) = [ q i 



(22) 



(23) 
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A compact way of expressing EqfJT) Eqf22] and Eqf23] is 
given by 

Finally we can write 

(c k (t)ct(t)) = R x (e k )(B k (t)Bl(t))RUe k ) 

= R x (0k)S z (-2A k t)(B k Bl)Si(-2A k t)Rl(d k ) 
= Rx(9k)S z (— 2A k t)R x (9' k — 9 k ){B' k B' k ) 
xRt(9 k -6 k )St(-2A k t)Rt(9 k ). (25) 

Notice that the above expression is the same if we consider 
S z (— 2A k t), with Afc the single particle spectrum of EqQ] or if 



we consider S z (— 2A k t), with A k the single particle spectrum 
of the Ising Hamiltonian without the DM term. This can be 
seen with a direct calculation. For example, one entry of the 
above correlation matrix is given by 

2<c t _ fc (t)c_ fc (t)) = 

Ex + E 2 + (E 2 - Ex) cos 9 k cos (9' k - 9 k ) 
+ (Ex - E 2 ) cos [t (A* + A_ fc )] sin0 fc sin (9' k - 9 k ) , (26) 

where Ex = |(1 + ^) and E 2 = |(1 - ^^). The argu- 
ment of S z appears only in the argument of the trigonometric 
function in such a way that the DM contribution is irrelevant 
(see Eq[7]i. This proves that the time evolution of the cor- 
relation matrix in Eqj25] with respect to the ground-state of 
EqQ] is the same as the time evolution of the correlation ma- 
trix with respect to the Ising Hamiltonian in transverse field, 
with respect to an excited state of the Ising Hamiltonian. 



